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We present a new model of large-scale multilayer convection in solar 
type stars. This model allows us to understand such self-similar structures 
observed at solar surface as granulation, supergranulation and giant cells. 
We study the slow-rotated hydrogen star without magnetic field with the 
Q spherically-symmetric convective zone. The photon's flux comes to the con- 

^ vective zone from the central thermonuclear zone of the star. The interaction 

of these photons with the fully ionized hydrogen plasma with T > 10 5 -fT is 
i — | carried out by the Tomson scattering of photon flux on protons and elec- 

trons. Under these conditions plasma is optically thick relative to the Tom- 
son scattering. This fact is the fundamental one for the multilayer convection 
formation. We find the stationary solution of the convective zone structure. 
This solution describes the convective layers responsible to the formation of 

2 the structures on the star's surface, 

-i— > 
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The systematic extreme ultraviolet and X-ray emission observations from 
Skylab station, Yohkoh, SoHO and Trace satellites give us the very inter- 
esting images of solar corona. After the previous images modifications (the 
partial gain of some interesting details) we can see large-scale corona struc- 
tures around the solar disk and these structures are not associated with 
. ^ active regions (Priest et al., 1990; Chertoc, 2002). The structures are similar 

to standard coronal loops that connected separate active regions together 
03 (Beck, 1998), but their "foots" lean on the photosphere out of active regions. 

These regular structures cover the hole solar disk as the more large-scale 
chromocpheric network. The lifetime of such loops is about a week for rel- 
atively small loops with length approximately equal to 2 • 10 4 km and with 
plasma concentration approximately equal to 10 15 m~ 3 but for the most great 
loops with length approximately equal to 3 • 10 5 km and with plasma con- 
centration approximately 2 • 10 14 m~ 3 the lifetime is about some months. 
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These large-scale structures (chains, loops) are observed for some years. We 
see that all the observed loop systems associated with quiet Sun permanent 
exist as a regular part of solar corona. It's necessary note that photosphere 
and chromosphere have regular structures such as grains, supergrains and 
giant grains. The giant grains are discovered by the helioseismology's meth- 
ods (Bee, 1998). These giant grains have the regular structure, their sizes 
are about 3 ■ 10 5 /cm with regular plasma speeds of lOOm/s. It's well known 
that granulation, chromospheric network, supergranulation and giant loops is 
the consequence of under-photosphere convective zone existence. The solar- 
like stars photospheres have the similar structures as Sun has: grains and 
supergrains. 

In this paper we present the simple model of the hydrogen star convection 
zone. 

The necessary condition of free convection (rises in plasma layers with 
thickness of only some times smaller then solar radius) is the Schwarzschild 
criterion - the specific entropy of plasma decreases with moving away from 
the star center. Such convection will develop when the temperature inside of 
small convective volume (convective cell) decreases slower than the temper- 
ature decreases in neighboring plasma (Rudiger, 1989). 

We use the Schwarzschild criterion later in our paper. 

If solar convection is laminar so such processes as granulation, chromo- 
sphere network and supergranulation may exist in the convective layers of 
different thickness. Therefore solar convective zone consists of the three lay- 
ers at least. 

Under large-scale laminar convection conditions the small-scale turbulent 
convection appears owing to development of different plasma instabilities. 
The ejections of matter at granulation and supergranulation scales are con- 
nected more than likely with instabilities dynamics. The lifetime of these 
ejections is small and they can't change the regular structure of quiet Sun 
but these ejections outline this structure by effective way. 

In this paper we propose the laminar convection model. In this model all 
the structures - granulation, chromosphere network, supergranulation and 
regular large-scale coronal structure - are examined as result of laminar con- 
vection action. We suppose that these structures are the different realizations 
of solution spectrum. 

Let's consider the main assumptions of this model: The convective zone 
is the layer with the spherically symmetry distribution of plasma around the 
radiative transfer energy zone. In this layer the condition of real hydrostatic 
equilibrium is carried out. 

In this paper we consider a case when the layer under study consists of 
the ionized hydrogen plasma only (protons and electrons). This consideration 
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may not be applied for the Sun atmosphere conditions but it significantly sim- 
plifies the mathematical description of convection and allow us understand 
the mechanism of convection zone structures formation. In our model this 
layer is open system through which the energy flux moves upwards. So let's 
consider that the plasma conditions we can described as polytropic equation: 



N_ 
No 




(1) 



where n is the polytropic index, N and T are the plasma concentration 
and temperature. These values are Nq ~ 5 ■ 10 27 m~ 3 and To ~ 2 • 10 6 K at the 
bottom border of the layer. The layer thickness (the convective zone depth) 
is approximately 0.3R Q , where R Q is the solar radius. The energy emission 
come to convective zone, the temperature T r T near the bottom border 
of the convective zone. This flux is the reason of the development of laminar 
convection. 

The emission and ionized plasma interaction is carried out by photon 
scattering on electrons and protons in case where the photon energy don't 
exceed the value kT. The time of energy transmission from photons to plasma 
don't excess the value (Kaplan, Tsytovich, 1973): 

if kT r <C m p c 2 , where k is the Boltzmann's constant, a? is the Thomson 
probability section of scattering, m p is the electron mass, e r = ^-T r 4 , Ob is 
the Stephan - Boltzmann constant. If T r = To then to ~ 0.1s. 

The distance of free run for photons is equal to A~ (err-AT) -1 . If N — Nq 
then Asi3m. 

In A 3 volume plasma and emission are in thermodynamic equilibrium 
almost because of the radiation is connected with matter. 

The thermal conductivity mechanism is made available by the next pro- 
cesses in our case. The plasma (heated by radiation in value A 3 ) loses the 
energy by bremsstrahlung. The speed of these losses is e ep = 1.6- 10~ 40 iV 2 v / T 
J ■ m~ 3 s _1 . The characteristic time of this process is equal to (Kaplan, 
Tsytovich, 1973): t 2 = 2.6 • 10 17 VTN^s. If N = N ,T = T we derive 
t 2 ~ lCT 6 s. At the other hand the bremsstrahlung heats up the electrons in 
the vicinity of the volume A 3 . The time taken for this process 

3m e c 

where E = 3/2NkTm e /m p is the electron energy density. If T < T < T 
and N* < N < N , N* = 2.5 • 10 26 m~ 3 we find 20s < h < 160s. 
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The characteristic rate of thermal conductivity is equal to v x = ^. So 
the thermal conductivity coefficient for the process described (for the same 
order of magnitude) is equal to: 

X = v x \ = — (4), 

where A is the thickness of the shell warmed up. 

The convectional energy transfer is carried out thanks to macroscopic 
transports of the value A 3 . The temperature inside the volume A 3 is higher 
then the plasma temperature in the layers which are situated higher then the 
bottom border of the convective zone, see Section 1. Thus the Archimedean 
raising force acts on this value and gives him the acceleration where 
g is the free fall acceleration on the bottom border of the convective zone, 
AT = T -T,T < T . 

The flotation process is retarded by viscosity. In our case the viscosity 
is the consequence of the Tomson scattering. The value A 3 is full of plasma 
and radiation. When this value moves the radiation is scattered by electrons 
of neighboring plasma. Thanks to the scattering the equalization of elec- 
tron momentum takes place inside the volume A 3 and outside of one. This 
viscosity they called radiation viscosity. It characterized by the viscosity 
coefficient 

V= l^N (5) - 

If N = N Q then v rs 6 • 10 9 m 2 /s. This value is similar to the value esti- 
mation taken from the analysis of observations. The floating is ended when 
the raising force is in equilibrium with viscosity forces. The characteristic 
time of convective floating is equal to 

to = 



9\W 

where A is correspond to characteristic scale of the convective layer (the 
mixing length). If A ~ 2 ■ 10 8 m, T = To and ^ ~ 1, g — 2g Q , where 
g Q ~ 274m/ s 2 is the gravity force acceleration on the solar surface then we 
have t 2 ~ 0.05s. So at the bottom border of the convective zone the relation 
t\ > t 2 is taken place. In this case the convective transfer is more effective 
then the heat conduction. 

Near the top border of the convective zone N — 4 • 10 22 m~ 3 and A* = 
(a^N)" 1 w 10 3 km. In this case we can ignore the Tomson effect. The 
radiation of plasma propagates free up to solar photosphere. 
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In the section 1. we give the solutions of stationary convective zone 
structures in the hydrodynamics approximation with the heat conduction 
(4) and viscous (5) coefficients. 

These solutions have the solitary wave structure and describe the model 
of multi-layer convection. All the convective cells have the torus contour. 

1 THE EQUATIONS OF THE STATIONARY CONVECTIVE ZONE 
STRUCTURE 

The set of simultaneous equations for the spherically symmetric station- 
ary convective zone which rotates about z axis (because of the hydrodynamics 
approximation is correct) have the form: 

(v,V)v= V ^ + ^) + g + u . Vv-2[v,Cj) (6) 
p 

is the motion equation, where cJ is the angular velocity of convective zone 
rotation, p is the plasma density, p is the plasma pressure, p r is the pressure 
of radiation. 



(v,V)T = X -AT (7) 
is the heat conduction equation, 

V(p + p r )+ P g = (8) 
is the hydrostatics equilibrium equation, 

dp r (JtN 1 



dr c Anr 2 

is the radiation transfer equation outside of the volume A 3 , 



(9) 



= 47rr 2 e ep (10) 
is the bremsstrahlung of plasma equation inside the volume A 3 , 

dM = 4irr 2 p-dr (11) 
is the mass conservation equation. 

In the equation (8) we don't take the density of radiation p r because of 
p r <C p in solar-like stars. 

The set of simultaneous equations (8-10) have used by A.S.Eddington in 
1926 (Eddington, 1926). 
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Let p = NkT,N = N(r),T = T(r) and state of plasma is described by 
the polytropic equation (1). Let's take the variable x — where 

C= [ -r-FT ~ ) ~ 13 • W 5 km > R Q . 
KAirGrripPo/ 

Then function r = ^ with equations (8) - (11) can be transformed to the 
following equation 

(n+l)- 2 {x 2 T x ) x = -T n + ar 2n+ K (12) 
where index x means the differentiation with respect to x and 

a = aT£ep(iVo ' To) « 10 17 . 
AirGiripPo 

Let velocity vector v has the {V, W 7 , Z} components in spherical coordi- 
nate system. The vector of angular velocity u has the following components: 
{ujcosO, —usin9, 0, 0}. 

Let <C 1. From the equations of the structure 

(v, V)v = v • Av — 2[v, uj] 

(v,V)T = x -AT, (13) 
one can find the equation for V(x) and t(x): 

v 1 



VV X = ---- (x 2 V x 
C x 2 



Vr. = l±- (x 2 ^)* (14) 



The equations (12) and (14) are simplified when we assume that the compo- 
nent of velocity V is decreased with the depth. This condition is in agreement 
with solar observations: the plasma spread out velocity in the photosphere 
decreases with the scale increasing from grains to giant grains. 
We choose the solution in the next form: 

V=- ( v (15) 

where a is the free parameter. 

This permits us to simplify the equation (14) and transform it to the 
following: 
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(n + 1)0t x = -r 2n+1 (l - ar n+1 / 2 ), (16) 

where f3 = 7.5, a = ^a. The equation (16) has different solutions for the dif- 
ferent values of n. Let's choose the value n (use the Schwarzschild criterion). 
According to this criterion the temperature inside the small element A 3 has 
to decrease with increasing of the distance from the star center slower then 
decreasing of plasma temperature occurs. The plasma is in the hydrostatic 
equilibrium and the radiation is absence. 

Substitute p r = 0,p = NkT and p = m p ■ N to (8). 

Use the polytropic equation (1), we find the relative change of the plasma 
temperature (radiation doesn't take into account) 

,Tx, = m p gC _ x 
't |0 (n + l)kT T ' 

In the volume A 3 the temperature changes according to (16). Also let's 
take into account a»l and r < 1. Then ar n+1 / 2 3> 1 and relative change 
of the temperature inside the volume A 3 is approximately equal to: 

~ 01 r 3n+l/2 

'T 1 ~ (n+l)/3 

In our case of the evolution of the convective instability |^| > |^| the 
number f3 is evaluate as: f3 > ^T 3n+3 ^ 2 - In this case we have r n 3> (cty/r) -1 . 
Therefore (3 > (223a 2 )- 1 . 

At the other hand one can integrate the equation (16) because of ignoring 
the first member of the right part of the equation. 

Thus we obtain the next algebraic equation: 

r -3n-i/2_ 1 = «3n + l/2 

P 71+1 1 

Using this equation and the consideration that r™ 3> (ay/r)" 1 we find that 

P < 3 n"+/ 2 ^^- ^ ne va l ue °f polytropic index n (1) is necessary to choose 
as to make up the next unequality: 

_J_ 3n + l/2x-xo 

223a 2 1 n + 1 a 2 r { ' 

For the solar-like star we have ^ < r < 1 and 10~ 4 < x — x < j^. In this 
case the unequality (17) is realized for all n having the positive values. Let's 
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choose n = 3/2. Then the accurate solution of the equation (16) for t(x) one 
can find from the next algebraic equation: 



2 . . 1 1 . 1. a 3 / 2 / a 1/2 + l , « 1/2 r+lx 

^ - xo) = 3 (1 " ^ + a(1 - r } - T" r ^v^i - 

(18) 

We have taken into account that t(x ) = 1 here. 

The solution for r(x) has the solitary wave form. It's clear from the form 
of the equation (16). 

If « x / 2 r >> 1 we have the asymptotic solution 

For x — > x Q one can find that 

t « e -^ { """ o) ->■ 1. 

At last we can find the expression for the speed components W and Z. 
Then we examine the most simple case of the symmetric spreading out on 
the sphere surface when W = Z. Let's consider also that the angular velocity 
w we can take from the equation 

dW v 1 d ( ^dW\ 

V ™ +2W-C0- cos6 = -.-.£- (x 2( ^-) . (20) 
ox C x ox ^ ox / 

As follows from the equation (20) the convective zone rotates differently. 
Thanks to the convection the redistribution of the rotatory moment inside 
the star takes place. This effect is accurately studied in (Rudiger, 1989). 

Under conditions selected in our paper we can find the equation for W 
from the first equation of the set of simultaneous equations (14). His form 
becomes simple enough 

dW 2 _ v 1 dW 2 
~dT ~dP' 

if we change the e and 4> angular variables to / variable and dl = a/ 'de 2 + sin 2 edip 2 . 
Among the multiple numbers of solutions of the equation (21) there is peri- 
odic solution. This periodic solution has the next form: 

W = W -tg(W ^-x-(l-l )), (22) 

where W is the speed peak value W, the point l is situated at the radius 
x and is the start reading for I coordinate. On the surface of sphere with 



(19) 



(21) 
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radius x* plasma spreads out from l point. So our model is symmetric 
there are many points Iq^ on the surface of the sphere of radius x*. The 
distance between the neighboring points is equal to 2£ = x (lQ ti — lo,i-i) = 
• A/ = ^^2- Between these points there are two opposing plasma streams 
with velocities of opposite direction. These streams compensate each other 
at the distance £ from the each points. 

So all the surface of the radius x* breaks-down to the cells with diameters 
which are equal to £. All the number of these cells L we can calculate when 
we the surface square nx\ divide by the cell square 7r(£/2) 2 : L — 4 • (x*/l;) 2 . 
Then the velocity amplitude is equal to Wq = irv /2£ 2 x*\fL. The kinetic 
energy density e is proportional to Wq. So the convective streams have the 
spectral energy distribution e ~ L ~ e~ 2 . The solutions of the convective 
zone structures (21) and (22) describes the stationary convection when all 
zone of the convective energy transfer consists of the layers with the different 
thickness. Every convective cell have the torus form. These solutions of 
this important problem are made for the first time. From the equation (19) 
follows the next conclusion: the convective zone differently rotates. Thanks 
to the convection the rotation moment redistribution inside the star is taken 
place. This effect is studied in detail in (Rudiger, 1989). 

2 SUMMARUY AND CONCLUSIONS 

This model qualitatively describes the deep convective layers of the star 
under the supergrains layer. In case of the star's convective transfer it's 
important that plasma at these layers is the fully ionized. We don't study 
star's plasma at the highest convective under-photospheric layer where the 
turbulent processes are possible. In this turbulent layer there are necessary 
conditions for the generation of the long-scale magnetic field of the star. At 
the layers under this turbulent under-photospheric layer the convection is the 
stationary convection. 

Let's use the asymptotic solution (19) for the convective zone analysis. 
We have the convective zone consists of some layers with thickness of Aj, i — 
0, 1, 2, ... The temperature on the lower part of the layer's border is equal to 
Tj, on the top part is Tj». If we take the dependence of parameters from (15) 
and (19) on To into account so we can find the relation between the velocity 
and temperature at the bottom and top borders of the neighboring layers: 

Vi-JVi = (W^m/T^ 372 ^^) 5 (23) 

For the qualitative estimation let's substitude the characteristics of the 
convective layers associated with giant cells and supergrains into (23): 

Vb = 10m/s,T ~ 2 • 10 6 K, A ~ 3 • 10 5 A;m 
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Figure 1: The convective zone structure 



Vi = lOOm/s, Ti ~ 10 6 fT, Ai ~ 3 • 10 4 km. 

In this case we obtain that To* ~ 0.4Ti* and the temperature on the top 
border of the layer Ao is smaller than the temperature on the top border of 
the layer Ai < A . So we can see that Ai torus are situated into A torus. 

This qualitative analysis of the formulae (23) allows us to make the con- 
clusion about relatively disposed convective layers in the hydrogen star. The 
layers are put one into another as we can see at the Figure 1. 
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In (Rozgacheva et al., 2003; Rozgacheva et al., 2004 ) was shown, that 
the torus typical scales may form the geometric progression. This fact is 
one of our model tests. Such geometric progression that described stationary 
structures at Solar surface is probably form the fractal set. 
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